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parallelized  C- language  version  of  the  code  has  been  developed  and  tested  using  a  distributed 
memory  transputer  array  consisting  of  21  T-800/20  MHz  processors. 

The  FDTD  code  has  been  used  to  calculate  currents  induced  in  the  human  body  and  the 
specific  absorptions  for  the  various  organs  for  exposure  to  vertically  polarized  electromagnetic 
pulse  (EMP)  using  both  the  5628-  and  45,024-cell  representation  of  the  human  body.  Also  the 
code  has  been  modified  to  allow  a  linearly  expanding  grid  representation  of  the  biological 
bodies.  Allowing  a  high-resolution  representation  of  the  body  surface  including  the  eyes  and 
the  gonads,  this  latter  improvement  would  be  most  helpful  in  extending  the  SAR  calculations 
to  microwave  frequencies. 
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AUTOMATED  RADIOFREQUENCY  RADIATION  DOSIMETRY 


INTRODUCTION 

The  interaction  of  radiofrequency  (RF)  and  microwave  (MW)  electromagnetic  (EM)  radiation 
with  biological  tissues  is  of  increasing  importance  from  the  standpoint  of  health  and  safety.  There 
is  considerable  literature  devoted  to  the  study  of  RF  and  MW  bioeffects  based  on  both  animal 
experimentation  and  theoretical  models.  Typically  bioeffects  are  correlated  with  mass-normalized 
rates  of  RF  or  MW  energy  absorption  (specific  absorption  rates  (SARs))  expressed  in  watts  per 
kilogram  (W/kg). 

The  American  National  Standards  Institute  (ANSI)  has  issued  a  guidance  standard  1 1  ]  for 
radiofrequency  radiation  (RFR)  exposure  level  needed  to  cause  behavioral  effects  in  laboratory 
animals.  This  standard  is  based  on  a  threshold  whole- body-averaged  SAR  of  4  W/kg.  The  ANSI 
guidance  standard  incorporates  a  safety  factor  of  10,  hence  whole-body-averaged  SAR  levels 
below  0.4  W/kg  averaged  over  a  6-min  period  are  considered  safe.  Recognizing  the  highly 
nonuniform  nature  of  SAR  distribution,  including  some  regions  where  there  may  be  fairly  high 
local  SARs,  the  ANSI  standard  further  stipulates  that  local  SAR  in  any  1  g  of  tissue  must  not 
exceed  8  W/kg. 

The  main  focus  of  our  project  has  been  to  develop  anatomically  based  EM  models  of  the 
human  body  to  obtain  quantitative  information  on  SARs  for  the  various  parts  of  the  body  for  expo¬ 
sure  to  RFR  for  a  variety  of  exposure  conditions.  Another  thrust  of  this  project  was  to  adapt  the 
numerical  SAR  algorithms  to  computing  workstations  and  to  distributed  memory  parallel  proces¬ 
sors  that  are  becoming  more  affordable  and  hence  more  readily  available  to  many  potential  users  of 
these  codes. 


FINITE-DIFFERENCE  TIM* -DOMAIN  METHOD 

The  most  promising  of  the  methods  for  calculation  of  SAR  distributions  at  the  present  time  is 
the  finite-difference  time-domain  (FDTD)  method  [2-8].  This  method  was  first  proposed  by  Yee 
[9]  and  later  developed  by  Taflove  and  colleagues  [10-12],  Umashankar  and  Taflove  [13],  Holland 
[14],  and  Kunz  and  Lee  [15].  Recently  it  has  been  extended  for  calculations  of  the  distribution  of 
electromagnetic  fields  (EMF)  in  a  man  model  for  incident  plane  waves  [2-6]  and  for  exposures  in 
the  near  fields  [7,  8],  In  this  method,  the  time-dependent  Maxwell's  curl  equations 


V  x  E  =  -p  ~~  ,  VxH=oE+£^ 

dt  dt 


are  implemented  for  a  lattice  of  cubic  cells.  The  components  of  E  and  H  are  positioned  about  a 
unit  cell  of  the  lattice  (Figure  1)  and  evaluated  alternately  with  half-time  steps.  The  goal  of  the 
method  is  to  model  the  propagation  of  an  EM  wave  into  a  volume  of  space  containing  dielectric 
and/or  conducting  structures  by  time  stepping,  i.e.,  repeatedly  implementing  a  finite-difference 
analog  of  the  curl  equations  at  each  cell  of  the  corresponding  space  lattice.  The  incident  wave  is 
tracked  as  it  first  propagates  to  the  dielectric  structure  and  interacts  with  it  via  surface  current  exci¬ 
tation,  spreading,  penetration,  and  diffraction.  Wave  tracking  is  completed  when  a  sinusoidal 


Figure  1 .  A  unit  cell  of  the  Yee  lattice  indicating  the  positions  for 
the  various  field  components. 


steady-state  behavior  for  E  and  H  fields  is  observed  for  each  lattice  cell.  A  second  approach  is  to 
illuminate  the  body  with  a  pulse  of  radiation  in  the  time  domain.  If  the  medium  is  nondispersive 
(i.e.,  its  dielectric  properties  do  not  change  with  frequency),  this  method  has  the  advantage  of  giv¬ 
ing  the  performance  of  the  medium  at  various  frequencies  from  the  inverse  Fourier  transform  of  the 
calculated  fields.  Though  Spiegel  and  his  collaborators  [2-3]  have  followed  this  approach,  we 
have  used  the  sinusoidal  variation  of  the  illuminating  fields  since  the  dielectric  properties  of  the  tis¬ 
sues  are  highly  dependent  on  frequency. 

The  Maxwell's  equations  (1)  can  be  rewritten  as  difference  equations.  We  assume  the  media 
to  be  nonmagnetic,  i.e.,  (I  =  p.Q  and  a  grid  point  of  the  space  defined  as  (i,  j,  k)  with  coordinates 
(i8,  j8,  k5)  where  8  is  the  cell  size. 

Upon  substituting 


E  = 


E 

2 


(2) 


the  equations  in  difference  form  for  Ez  and  Hz,  for  examples,  can  be  written  as: 


_n+l  _n 

Ez  (i,  j,  k  +  1/2)  =  CAZ  (i,  j,  k  +  1/2)  Ez  (i,  j,  k  +  1/2)  + 


2 


(3) 


CBz(i,  j,  k  +  1/2)  • 


H^+1/2  (i  +  1/2,  j,  k  +  1/2)  -  Hy+V2  (i  -  1/2,  j,  k  +  1/2) 

+  H^+1/2  (i,  j  -  1/2,  k  +  1/2)  -  H;+1/2  (».  j  +  1/2,  k  +  1/2) 


Hz+,/2  (i  +  1/2,  j  +  1/2,  k)  =  h"  1/2  (i  +  1/2,  j  +  1/2,  k) 

E"  (i  +  1/2,  j  +  1,  k)  -  e"  (i+  1/2,  j,  k) 

+ 

+  Ey  (i,  j  +  1/2,  k)  -  Ey  (i  +  1,  j  +  1/2,  k) 


where 


CAz(i,j,k+  1/2)  = 


gz  (i,  j,  k  +  1/2)  5t 
2a,  (i,  j,  k  +  1/2) 


_-l 

az  (i,  j,  k  +  1/2)  8t 
1  +  2e^  (i,  j,  k  +  1/2) 


(5) 


CBZ  (i,  j,  k+  1/2)=- 


e  (i,  j,  k  +  1/2)  + 


a  (i,  j,  k  +  1/2)  5t 


-l 


(6) 


where  the  superscript  n  denotes  the  time  n  8t  in  terms  of  the  incremental  time  or  time  step  8t.  The 
time  step  8t  is  determined  by  the  cell  size  and  must  satisfy  the  stability  condition  [10] 


vmax  ^ 


5\z 


-V— 

8y2  8z2 


-,-1/2 


(7) 


where  v  x  is  the  maximum  wave  phase  velocity  within  the  model.  For  the  cubic  cell  model  8x  = 
8v  =  8z  =8,  the  following  relationship  is  usually  taken 


8t  = 


2v, 


max 


(8) 


where  v  is  the  maximum  velocity  of  EM  waves  in  the  interaction  space.  We  have  taken  vnm 
corresponding  to  velocity  of  EM  waves  in  air.  The  e  and  o  may  be  different  for  Ex,  Ev.  and  Ez  in 


3 


the  inhomogeneous  media,  because,  as  shown  in  Figure  1,  the  components  of  E  are  not  positioned 
at  the  same  point  for  the  same  cell  in  Yee's  algorithm.  Thus  we  have  e*,  Ev,  e  ,  o  ,  o  ,  and  o  for 
each  of  the  cells. 

A  basic  problem  with  any  finite-difference  solution  of  Maxwell's  equations  is  the  treatment  of 
the  field  components  at  the  lattice  truncation.  Because  of  the  limited  computer  storage,  the  lattice 
must  be  restricted  in  size.  Proper  truncation  of  the  lattice  requires  that  any  outgoing  wave  disap¬ 
pear  at  the  lattice  boundary  without  reflection  during  the  continuous  time  stepping  of  the  algorithm. 
An  absorption  boundary  condition  for  each  Field  component  is  therefore  needed  at  the  edge  of  the 
lattice.  In  our  formulation,  the  absorption  boundary  conditions  in  second  approximation  derived 
by  Mur  [  16]  are  used.  The  condition  for  Ez  at  i  =  1,  for  example,  can  be  written  as: 


_n+l  _n-l 

Ey  (l.j,  k+  1/2)  =  -E,  (2,j,  k+1/2) 


^-.-i  bt  -  §  n+1  _n-l 

+  c^srrs  l E/  (2’j-k  + 1/2)  +  Ez  nj,k  +  i/2) 


( 


yc,  5t  +  S 


(C,  8t) 

5(C.  8t  +  5) 

/ 


e"  (l,j,k+  l/2)  +  E/n(2,j,k+  1/2)) 


<C  fitf 
25<(' .  6t  +  oi 


E."(l.j  +  1,  k  +  1/2)  +  Ez  (1,  j  -  1,  k  +  1/2)  +  Ez  (2,  j  +  1,  k  +  1/2) 
+  Fvn(2.j-  1,  k  +  1/2)  +  Ez  (1,  j,  k  +  3/2)  +  Ez  (1,  j,  k  -  1/2) 

+  e"  <2.  j.  k  +  3/2)  +  Ez  (2,  j,  k  -  1/2) 


(9) 


where  C  _  is  the  phase  velocity  of  the  medium  at  the  lattice  truncation. 

As  a  pan  of  this  project,  several  steps  have  been  taken  to  improve  the  efficiency  of  the  FDTD 
code.  Since  a  major  problem  with  this  code  is  the  large  amount  of  memory  storage  that  is  needed, 
we  have  reduced  the  parameters  that  need  to  be  stored.  First,  a  ,  a  <5  and  e  ,  e  ,  and  e  that 
were  previously  stored  for  all  N  cells  at  half-cell  intervals  have  been  replaced  by  single  volume- 
averaged  values  of  a  and  e  for  each  of  the  cells,  resulting  in  considerable  reduction  iiuhe  number 
of  parameters  that  need  to  be  stored.  This  step  has  also  resulted  in  the  intermediate  propertv- 
dependent  quantities  such  as  CAX,  CB^,  CA  ,  CB  ,  CA  ,  and  CB  being  replaced  bv  CA  and  CB 
tor  each  of  the  cells  -  leading  once  again  to  a  considerable  reduction  in  memorv.  The  steps  taken 
to  date  have  reduced  the  memory  storage  requirement  by  roughlv  one  half  and  computation  time  for 
a  typical  run  for  N  =  100,000  cells  using  CRAY  XMP  from  8  to  5  min 


4 


With  these  memory  and  time-saving  steps  taken,  the  test  runs  for  currents  induced  in  the 
human  body  for  exposure  to  electromagneuc  pulse  (EMP)  have  shown  negligible  difference  in  the 
results. 


AN  ANATOMICALLY  BASED  INHOMOGENEOUS  MODEL  OF  THE  HUMAN  BODY 

As  previously  described  in  Sullivan  et  al.  [4,  5],  the  model  for  the  human  body  is  developed 
from  information  in  the  book,  A  Cross-Section  Anatomy ,  by  Eycleshymer  and  Sehoemaker  117]. 
This  book  contains  cross-sectional  diagrams  of  the  human  body  which  were  obtained  by  making 
cross-sectional  cuts  at  spacing  of  2.54  cm  (1  in.)  in  human  cadavers.  The  process  for  creating  the 
data  base  of  the  man  model  was  the  following:  first  of  all,  a  0.635-cm  (1/4  in.)  grid  was  taken  for 
each  single  cross-sectional  diagram  and  each  cell  on  the  grid  was  assigned  a  number  corresponding 
to  1  of  the  14  tissue  types  given  in  Table  1,  or  air.  Thus  the  data  associated  with  a  particular  layer 
consisted  of  three  numbers  for  each  square  cell:  x  and  y  positions  relative  to  some  anatomical  ref¬ 
erence  point  in  this  layer,  usually  the  center  of  the  spinal  cord,  and  an  integer  indicating  which  tis¬ 
sue  that  cell  contained.  Since  the  cross-sectional  diagrams  available  in  Eycleshymer  and 
Sehoemaker  [17]  are  for  somewhat  variable  separations,  typically  2. 3-2. 7  cm,  a  new  set  of  equi- 
spaced  layers  were  defined  at  0.635-cm  (1/4  in.)  intervals  by  interpolating  the  data  onto  these  lay¬ 
ers.  Because  the  cell  size  of  0.635  cm  results  in  a  model  with  a  large  number  of  cells  that  is  diffi¬ 
cult  to  accommodate  in  the  memory  space  of  commonly  available  computers,  the  proportion  of 
each  tissue  type  was  calculated  next  for  somewhat  larger  cells  of  size  1.27  cm  (1/2  in.)  or  2.54  cm 
(1  in.)  combining  the  data  for2x2x2  =  8or4x4x4  =  64  cells  of  the  smaller  dimension. 
Without  changes  in  the  anatomy,  this  process  allows  some  variability  in  the  height  and  weight  of 
the  body.  We  have  taken  the  final  cell  size  of  1.31  or  2.62  cm  (rather  than  1/2  in.  or  1  in.)  to 
obtain  the  whole-body  weight  of  69.6  kg  for  the  model.  The  number  of  cells  either  totally  or  par¬ 
tially  within  the  human  body  for  the  two  models  are  45,024  and  5,628,  respectively. 


CURRENTS  INDUCED  IN  THE  HUMAN  MODEL  FOR 
EXPOSURE  TO  ELECTROMAGNETIC  PULSE 

Though  not  assigned  as  a  task  on  this  project,  we  have  used  the  FDTD  method  to  calculate 
the  internal  fields  and  the  induced  current  densities  in  anatomically  based  models  of  the  human 
body  with  5,628  or  45,024  cubical  cells  of  dimensions  2.62  and  1.31  cm,  respectively,  for  expo¬ 
sure  to  EMP.  High-intensity  electric  and  magnetic  fields  associated  with  an  EMP  have  resulted  in  a 
need  to  understand  its  coupling  to  the  human  body  to  assess  potential  health  effects.  Electric  fields 
as  high  as  50-100  kV/m  and  magnetic  fields  that  are  approximately  1/377  times  lower  have  been 
measured  close  to  EMP  generators.  The  fields  are,  of  course,  considerably  lower  at  larger  dis¬ 
tances  from  these  generators.  The  EMPs  are  marked  by  rise  times  on  the  order  of  10-30  ns  with 
pulse  durations  on  the  order  of  100-300  ns. 

Approximating  the  human  body  by  3  cylinders  of  uniform  conductivity.  1  vertical  and  2  hori¬ 
zontal  to  represent  shoulders,  Grrmhaug  has  calculated  currents  induced  for  exposure  to  an  EMP 
1 18|.  Guy,  on  the  other  hand,  .<as  used  a  model  consisting  of  12  cylindrical  sections  each  repre¬ 
senting  approximately  15-cm  (5.9  in.)  height  of  the  human  body  1 19],  These  sections  are  loaded 
w  ith  parallel  resistance-capacitance  (RC)  circuits  representing  the  measured  impedance  for  a  human 
subject  for  the  corresponding  sections  of  the  body.  Approximating  the  incident  EMP  hv  an 
expression  Eft)  =  K[em  -  e^],  currents  for  various  frequency  components  (75  frequencies 
covering  the  range  100  Hz  to  300  MHz)  are  calculated  using  the  NEC  2D  code  detailed  in  Burke 
and  Poggio  [2()|.  The  induced  current  in  the  time  domain  for  the  minimum  area  section  rep¬ 
resenting  the  feet  and  ankles  is  then  obtained  from  the  Fourier  transform  of  the  currents  calculated 
for  the  various  frequencies. 
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TABLE  1.  TISSUE  PROPERTIES  USED  AT  SOME  REPRESENTATIVE  FREQUENCIES 


Tissue  Type 

Mass  density 
1,000  kg/m3 

27.12  MHz 

°S/m  ^ 

100  MHz 

aS/m 

350  MHz 

°S/m 

Air 

0.0012 

0.0 

1.0 

0.0 

1.0 

0.0 

1.0 

Muscle 

1.05 

0.74 

106 

1.0 

74 

1.33 

53.0 

Fat,  bones 

1.20 

0.04 

29 

0.07 

7.5 

0.072 

5.7 

Blood 

1.00 

0.28 

102 

1.1 

74 

1.2 

65.0 

Intestines 

1.00 

0.29 

60 

0.55 

36 

0.66 

26.5 

Cartilage 

1.00 

0.04 

29 

0.07 

7.5 

0.072 

5.7 

Liver 

1.03 

0.51 

132 

0.62 

77 

0.82 

50.0 

Kidney 

1.02 

0.79 

209 

1.0 

90 

1.16 

53.0 

Pancreas 

1.03 

0.69 

206 

1.0 

90 

1.16 

53.0 

Spleen 

1.03 

0.69 

206 

0.82 

100 

0.9 

90.0 

Lung 

0.33 

0.17 

34 

0.34 

75 

1.1 

35.0 

Heart 

1.03 

0.64 

210 

0.75 

76 

1.0 

56.0 

Nerve,  brain 

1.05 

0.45 

155 

0.53 

52 

0.65 

60.0 

Skin 

1.00 

0.74 

106 

0.55 

25 

0.44 

17.6 

Eye 

1.00 

0.45 

155 

1.9 

85 

1.9 

80.0 

Since  these  models  are  fairly  gross  approximations,  we  have  used  the  anatomically  based 
model  of  a  human  to  calculate  currents  induced  in  the  human  body  for  exposure  to  a  few  represen¬ 
tative  EMPs.  The  FDTD  technique  [5-7]  was  used  to  calculate  the  internal  fields  (E)  and  the  cur 
rent  densities  (J)  for  each  of  the  cubical  cells  of  dimensions  2.62  cm  (5,628  cells)  or  1.31  cm 
(45,024  cells)  that  are  used  to  describe  the  human  body.  As  previously  described  [5-7],  the  vol¬ 
ume-averaged  dielectric  properties  and  mass  densities  were  obtained  from  the  contents  of  these 
subvolumes  in  terms  of  the  16  tissues  and  air.  The  information  on  the  tissue  contents  were 
obtained  with  a  resolution  of  0.635  cm  (1/4  in.)  from  the  anatomical  diagrams  [17].  The  models 
were  assumed  with  a  layer  of  dielectric  constant  e^.  =  4.2  of  thickness  2.62  cm  under  the  feet  to 
simulate  a  human  wearing  rubber-soled  shoes.  This  assumed  thickness  corresponds  to  the  width 
of  1  or  2  cells  for  the  2  models,  respectively. 

Being  a  time-domain  technique,  the  FDTD  method  is  ideally  suited  for  calculations  involving 
EMP.  Exact  pulse  shapes  can  be  modeled  by  using  small  time  increments.  This  method  does, 
however,  require  that  the  dielectric  properties  of  the  body  be  nondispersive,  i.e.,  independent  of 
frequency  in  the  band  representing  the  frequency  content  of  the  EMP.  This  requirement,  unfortu¬ 
nately,  is  not  met  for  the  human  tissues  which  are  known  to  be  dispersive.  A  rigorous  approach  to 
this  problem  using  the  FDTD  would  therefore  require  the  solution  of  Maxwell's  equations  coupled 
to  the  polarization  equations  in  the  time  domain  for  the  various  tissues.  Since  the  scope  of  the 
present  project  did  not  call  for  a  full-blown  effort  in  this  regard,  we  have  used  an  approximate 
model  that  uses  the  dielectric  properties  at  40  MHz.  This  approach  is  somewhat  justified  since  the 
predominant  component  of  the  induced  current  are  at  frequencies  close  to  40-50  MHz,  which  is 
near  the  resonant  frequency  of  an  adult  human  standing  on  a  ground  plane  [21], 

The  total  induced  currents  for  the  various  sections  of  the  body  and  the  specific  absorptions 
for  several  organs  are  given  by  Chen  and  Gandhi  [31]  for  two  representative  EMPs.  The  calcu¬ 
lated  results  for  the  induced  currents  are  compared  with  the  data  measured  for  a  human  subject 
(measured  data  courtesy  of  Mr.  Wayne  Hammer,  Naval  Surface  Warfare  Center,  White  Oak, 
Maryland).  The  calculated  results  are  in  excellent  agreement  with  the  measured  data  albeit  for  one 
human  subject  to  date. 

Using  the  graphics  available  with  the  IBM  model  520  workstation,  the  staff  at  IBM,  Salt 
Lake  City,  has  helped  us  create  a  movie  of  EMP  interaction  with  a  nominal  1-in.  model  of  the 
human  body  (actual  cell  size  =  2.62  cm).  Since  software  was  not  available  to  transfer  the  computer 
movie  onto  a  videotape,  pictures  of  sequential  frames  of  this  movie  were  used  for  our  presentation, 
"Current  Induced  in  an  Anatomically  Based  Model  of  a  Human  for  Exposure  to  Vertically 
Polarized  EMP,"  at  the  Twelfth  Annual  Meeting  of  the  Bioelectromagnetics  Society,  June  10-14, 
1990,  in  San  Antonio,  Texas. 

We  have  also  developed  a  simplified  procedure  to  calculate  the  coupling  of  the  EN‘P  to  the 
human  body.  In  this  procedure  the  response  of  the  human  model  is  calculated  for  a  narrow 
impi.'se  (width  of  0.218  ns)  of  RFR.  The  calculated  current  variations  for  the  various  sections  of 
the  body  (such  as  sections  through  the  ankles,  knees,  thighs,  bladder,  heart,  neck,  eyes,  and  oth¬ 
ers)  are  stored  on  a  computer  diskette.  By  convolving  the  calculated  impulse  response  with  a  pre¬ 
scribed  EMP  in  the  time  domain,  we  can  get  the  currents  for  these  same  sections  of  the  body  that 
would  have  been  induced  as  a  result  of  this  EMP.  This  procedure  has  been  tested  with  an  EMP 
provided  by  the  ALECS  facility,  Kirtland  Air  Force  Base,  New  Mexico  (Figure  4  [31  ]).  The  cur¬ 
rent  variations  obtained  for  the  various  sections  using  this  simplified  procedure  are  identical  to 
those  shown  in  Figure  6a-e  of  Chen's  and  Gandhi's  manuscript  [31],  the  latter  having  been 
derived  by  the  exact  method.  Since  the  convolution  integral  is  a  relatively  simple  step  that  can  be 
carried  out  with  a  personal  computer  (PC),  the  response  of  the  human  body  can  thus  be  obtained 
for  any  EMP  once  the  computer  run  for  the  impulse  response  has  been  made  and  data  supplied  on  a 
diskette. 
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Further  evaluation  of  the  previous  calculations  for  EMP  [31]  have  also  revealed  a  problem. 
Shown  in  Figure  2  is  the  Fourier  Spectrum  of  the  current  induced  for  the  section  through  the  knees 
taken  from  Figure  6b  of  Chen's  and  Gandhi's  report  {31].  Fairly  similar  Fourier  Spectra  have  also 
been  obtained  for  other  sections  of  the  body.  The  broad  bandwidth  of  the  current  components  (~ 
0-70  MHz)  casts  doubt  on  the  validity  of  having  assumed  frequency-independent  properties  e  and 
C  corresponding  only  to  40  MHz,  since  the  properties  over  the  pertinent  region  of  frequencies  can 
change  by  a  factor  of  2:1  or  more.  An  approach  to  resolve  this  problem  will  be  to  obtain  time- 
domain  formulations  of  the  tissue  properties  e  and  o  from  the  known  frequency-domain  variations 
and  use  these  to  calculate  the  impulse  response  which  can  then  be  convolved  with  the  prescribed 
EMPs.  This  approach  has,  however,  yet  to  be  tried  and  tested  for  accuracy. 


Figure  2.  Fourier  Spectrum  of  the  EMP-induced  current  for  the  section  through  the 
knees  (fairly  similar  Fourier  Spectra  have  also  been  obtained  for  currents 
for  the  sections  through  ankles,  thighs,  heart,  neck,  and  others). 


SPECIFIC  ABSORPTION  RATE  AND  INDUCED  CURRENT  DISTRIBUTIONS 
FOR  PLANE-WAVE  EXPOSURES  (100-915  MHz) 

We  had  previously  used  the  FDTD  method  to  calculate  local,  layer-averaged,  and  whole- 
body-averaged  SARs  and  internal  RF  currents  in  a  5,628-cell  (cell  size  =  2.62  cm)  anatomically 
based  model  of  the  human  body  for  plane-wave  exposure  20-100  MHz  [6]. 

Using  the  finer  resolution  (cell  size  =1.31  cm),  45,024-cell  representation  of  the  human 
body,  calculations  have  recently  been  extended  to  frequencies  of  120-200  MHz  in  10-MHz  steps 
and  frequencies  of  225,  250,  300,  350,  500,  700,  and  915  MHz.  Frequency-dependent  properties 
for  the  various  tissues  were  taken  from  the  literature  [22,  23].  This  work  is  described  in  detail  in 
the  repoit  by  Gandhi  et  al.  [32],  Special  attention  was  given  to  the  frequency  band  120-200  MHz 
to  see  if  we  could  predict  the  high  SARs  that  had  been  observed  for  the  neck  by  Stuchly  et  al.  [ 24] 
at  160  MHz  using  a  four-tissue  (skeleton,  brain,  lungs,  and  muscle)  heterogeneous  model  of  man. 
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In  the  work  by  Stuchly  etal.  [24]  only  three  frequencies  --  160,  350,  and  915  MHz  --  and  isolated 
human  model  had  been  studied.  By  taking  a  closer  frequency  interval  of  10  MHz  in  the  band  1 20- 
200  MHz,  we  wanted  to  pinpoint  the  frequencies  for  maximum  SARs  in  the  neck  both  for  isolated 
and  grounded  conditions  of  the  model.  From  our  calculations  for  the  isolated  human  model 
(irradiation  condition  used  by  Stuchly  et  al.  [24]),  the  highest  layer-averaged  SAR  is  found  for  the 
neck  section  not  at  160  MHr.  used  as  one  of  the  irradiation  frequencies  by  these  authors  (the  other 
two  irradiation  frequencies  were  350  and  915  MHz),  but  rather  at  200  MHz.  The  maximum  layer- 
averaged  SAR  for  the  neck  section  at  200  MHz  has  the  value  0.55  W/kg  for  an  incident  power 
density  of  1  mW/cm2.  Since  only  local  values  are  given  for  38  locations  in  the  body  in  the  article 
by  Stuchly  et  al.  [24],  data  was  not  available  for  comparison  with  the  layer-averaged  SARs  calcu¬ 
lated  at  160  MHz  or  for  any  of  the  other  15  frequencies  used  in  our  calculations.  The  highest  local 
SAR  calculated  for  the  center  of  the  neck  is  3.22  W/kg  at  200  MHz.  This  value  is  considerably 
higher  than  about  1.8  W/kg  that  was  measured  for  this  general  region  by  Stuchly  et  al.  [24]  at  160 
MHz  using  their  four-tissue  heterogeneous  model  (skeleton,  brain,  lungs,  and  muscle).  The  peak 
value  for  the  center  of  the  neck  in  our  16-tissue  anatomically  based  model  at  160  MHz  is  2.41 
W/kg.  This  value  compares  favorably  with  1.8  W/kg  measured  for  this  general  location  by 
Stuchly  et  al.  T24]  when  one  realizes  that  only  a  four-tissue  phantom  was  used  by  these  authors. 
While  no  reason  was  offered  by  Stuchly  et  al.  [24]  for  the  high  SARs  observed  in  the  neck  at  160 
MHz,  we  find  that  the  peak  RF  induced  currents  both  for  isolated  and  grounded  conditions  shift  to 
the  upper  half  of  the  body  for  frequencies  in  the  band  120-300  MHz  (see  Figure  9b- i  of  [32]) 
rather  than  the  lower  parts  of  the  body  and  the  legs  as  previously  reported  for  frequencies  smaller 
than  100  MHz  [6,  25].  Because  of  the  smaller  cross-sectional  area  of  the  neck,  highest  layer-aver¬ 
aged  and  local  SARs  are  therefore  calculated  for  the  neck  section. 

We  have  also  used  the  calculated  SAR  distributions  to  obtain  part-body  (head  and  neck, 
torso,  legs,  and  arms)  and  organ-averaged  SARs  (eyes,  brain,  heart,  lungs,  liver,  kidneys,  intes¬ 
tines).  Taken  from  the  report  by  Gandhi  et  al.  [32],  the  part-body  and  organ-averaged  SARs  are 
given  in  Figures  3  and  4,  respectively.  The  highest  part-body-averaged  SAR  for  the  head  and  neck 
region  for  the  isolated  condition  of  the  model  occurs  at  200  MHz  where  the  SAR  for  this  region  of 
the  body  is  182.5  mW/kg. 

In  Figures  4a  and  b  we  have  plotted  the  averaged  SARs  calculated  for  the  various  organs  in 
the  body  for  isolated  and  grounded  conditions  of  the  model,  respectively.  The  organs  considered 
are  brain,  eyes,  heart,  lungs,  liver,  kidneys,  and  intestines.  In  the  frequency  band  100-200  MHz 
the  highest  SARs  for  the  eyes  and  the  brain  occur  for  frequencies  that  coincide  with  the  highest 
absorption  frequencies  in  the  head  and  neck  region.  The  highest  absorption  frequencies  for  the 
eyes  and  the  brain  (as  well  as  the  head  and  neck  region)  are  200  and  150  MHz  for  isolated  and 
grounded  model  conditions,  respectively.  Because  of  the  increasingly  superficial  nature  of  EM 
deposition  at  higher  frequencies,  the  average  SARs  for  the  eyes  and  the  brain  also  rise  for  fre¬ 
quencies  above  350  MHz. 


MODIFICATION  OF  THE  FDTD  ALGORITHM  -  LINEARLY 
EXPANDING  GRID;  UNEQUAL  CELL  SIDES 

Recognizing  that  the  EM  deposition  at  microwave  frequencies  is  primarily  superficial,  e.g., 
in  the  eyes,  the  frontal  lobe  of  the  brain,  gonadal  area,  and  others,  we  have  modified  the  FDTD 
algorithm  to  use  a  linearly  expanding  grid  rather  than  a  uniform  grid  that  has  always  been  used  in 
the  past  for  this  method.  The  expanding  grid  algorithm  allows  cell-to-cell  expansion  factors  that 
can  be  different  along  the  three  axes.  Expansion  factors  a  on  the  order  of  1 .03  to  1 .08  have  been 
found  to  give  results  that  are  in  excellent  agreement  with  those  obtained  analytically  for  two- 
dimensional  (2-D)  homogeneous  and  layered  cylinders  and  for  three-dimensional  (3-D)  homoge¬ 
neous  spheres  [26].  Because  of  the  power  of  compounding,  the  cell  size  after  n-cells  is  6  an 
where  6mjn  is  the  minimum  cell  size  for  the  grid.  This  size  can  be  considerably  larger  thaiuhe 
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Figure  4.  Average  SARs  for  the  various  organs  of  the  body  for  an  incident  E-polarized 
power  density  of  1  mW/cm2.  Frequency  band  100-915  MHz. 


minimum  cell  size  that  is  used.  In  air  outside  the  body,  such  larger  cell  sizes  are  obviously  accept¬ 
able  and  are  still  much  smaller  than  X^IO  where  XQ  is  the  wavelength  in  air.  In  the  interior  of  the 
bodv  these  larger  cell  sizes  are  also  acceptable  since  the  SAR  is  fairly  negligible.  A  comparison  of 
the  fields  calculated  with  expanding  grid  algorithm  and  the  analytical  series  solution  [26]  for  a 
lossy  dielectric  sphere  (cr  =  46.  a  =  2.26  S/m)  of  diameter  21  cm  (8.4  in.)  at  3,000  MHz  along 
various  cuts  of  the  bodv  is  given  in  Figure  5.  The  agreement  can  be  seen  to  be  fairly  good.  Fairly 
excellent  agreements  have"  also  been" observed  for  this  sphere  at  other  frequencies  up  to  3,000 
MHz. 


W  e  have  also  modified  the  FDTD  algorithm  so  that  it  is  possible  to  use  different  cell  sizes 
along  the  three  axes.  In  future  studies  it  should  be  possible  to  use  parallelopiped  rather  than 
cubical  cells,  allowing  a  degree  of  flexibility  that  may  be  helpful  in  some  applications. 


ADAPTATION  OF  THE  FDTD  CODE  TO  COMPUTING  WORKSTATIONS 

A  major  task  on  the  project  was  to  adapt  the  FDTD  code  to  computing  workstations  that  are 
becoming  more  affordable  and  hence  more  readily  available  to  many  potential  users  of  these  codes. 
As  a  part  of  this  project  the  FDTD  code  has  been  adapted  and  used  with  the  following  worksta¬ 
tions:  Silicon  Graphics  IRIS  (memory  12  Mbytes),  APOLLO  RGB/DN  10,000  (memory  64 
Mbvtes).  and  IBM  320  (64  Mbytes).  In  addition  to  the  nominal  1-in  and  1/2-in  models  that  we 
have  used  for  the  results  described  in  the  earlier  sections,  a  nominal  1/4-in. -resolution  (actual  cell 
size  =  6.53  mm),  anatomically  based  model  of  the  human  body  has  been  recently  created  from  the 
data  base  that  was  previously  used  for  the  coarser  models.  This  model  is  nested  in  an  interaction 
space  of  2X5  x  1 10  v  62  =  1 .643.700  or  nearly  2  million  cells.  This  is  the  first  time  that  a  human 
body  model  with  such  a  high  resolution  has  been  run  anywhere.  For  preliminary  runs  we  have 
assumed  plane  wave  exposures  at  415.  1.500.  2.450.  and  3,0(X)  MHz  and  have  used  the  tissue 
properties  at  the  corresponding  frequencies.  APOLLO  10.000  Workstation  has  been  used  for  cal¬ 
culations  of  SAR  distributions.  Total  computation  limes  for  the  various  frequencies  are  as  follow’s: 
S  h  at  415  MHz:  5  ti.  20  min  at  1.500  MHz:  4  h.  40  min  at  2,450  MHz:  and  about  4  h  at  3,000 
MHz.  Because  of  the  need  to  cover  about  four  periods  of  the  incident  wave  to  obtain  convergence 
and  the  fact  that  the  time  step  8t  =  8, '2c  =  0.0104  ns  is  dependent  only  on  the  cell  size  and  not  on 
frequency,  fewer  iterations  arc  needed  at  higher  frequencies  --  hence  the  diminishing  computation 
:.me  at  these  frequencies 

Since  the  data  base  for  the  1/4-in.  model  is  less  than  satisfactory  having  been  obtained  by 
interpolation  from  anaiomic.il  sections  that  for  the  previous  models  were  available  only  for  separa¬ 
tions  on  the  order  of  2  3-2.7  cm.  it  is  obvious  that  a  better  model  is  needed  before  the  higher  fre¬ 
quency  results  are  taken  seriously.  Also  we  need  to  increase  the  separation  to  the  absorbing 
boundaries,  which  would  require  a  somewhat  larger  number  of  cells  for  the  interaction  space.  We 
ma\ .  of  course,  be  able  to  move  the  absorbing  boundary  behind  the  model  closer  to,  and  perhaps 
w  ithin.  the  body  since  the  fields  at  these  locations  are  negligible,  particularly  at  the  higher  frequen¬ 
cies. 


To  sum  up.  models  of  the  human  body  and  the  surrounding  space  (needed  for  truncation  of 
the  FDTD  code)  used  to  date  have  the  following  number  of  cells: 

84  x  38  x  26  =  82.992  cells  (cell  size  =  2.62  cm,  nominal  1") 

163  x  76  x  52  =  644, 176  cells  (cell  size  =  1.31  cm,  nominal  1/2") 

285  x  1 10  x  62  =  1,943,700  cells  (cell  size  =  0.655  cm,  nominal  1/4") 

The  memory  requirements  for  these  models  and  the  workstations  used  to  date  are  as  follows: 


0.0  L— 
-0.02 


0.08  0.13 

y(m) 


a.  yz  plane  at  the  center  (x  =  0) 


Lossy  dielectric  sphere 


i-iH 

r  L 

1.0  h  \ 
0.9  F-  . 


0.5  h 

0.4  h- 
I- 

0.3  t- 


0.1  - 


0.0  L— 

-0.02 


0.03  0.08  0.13  0.18 

v(m) 


b.  yz  plane  at  x  =  ^  =  2.625  cm 


f  igure  5.  Comparison  of  the  E-fields  calculated  from  linearly  expanding  grid  FDTD  algorithm 
with  the  analytic  results  (26|  for  the  various  cuts  of  a  lossy  muscle-equivalent  sphere  at 
3.000  MHz  (kQa  =  6.6).  Diameter  of  the  sphere  =  21  cm  (corresponding  to  the 
dimension  of  the  human  head).  Taken  from  reference  23,  dielectric  properties  of  the 
sphere  are  er  =  46,  c  =  2.26  S/m.  Expansion  factors  along  the  three  axes  are:  cxx  =  u/ 
=  1.08.  av  =  1.05. 
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1 .  82,992-cell  model  (less  than  12  Mbytes;  Silicon  Graphics,  IRIS,  APOLLO  10,000;  com¬ 
putation  time  approximately  7  s  per  iteration) 

2.  644,176-cell  model  (approximately  26  Mbytes;  APOLLO  10,000  and  IBM  520;  compu¬ 
tation  time  approximately  36  s  per  iteration  on  APOLLO  10,000  and  about  4  s  on  IBM 
520) 

3.  1,943,700-cell  model  (approximately  75  Mbytes;  APOLLO  10,000  with  extended  mem¬ 
ory;  computation  time  approximately  72  s  per  iteration) 


PARALLELIZATION  OF  THE  FDTD  CODE  AND  ITS  IMPLEMENTATION 
ON  A  DISTRIBUTED  MEMORY  TRANSPUTER 

s  work  done  by  Professor  Kris  Sikorski  and  his  student  Youn-Hsin  Young  is  described 
in  dc,.ui  [33].  As  a  result  of  this  task  the  FDTD  code  has  been  rewritten  in  a  parallelized  C-lan- 
guage  version  and  run  on  a  distributed  memory  transputer  array  consisting  of  21  T-800/20  MHz 
processors  interconnected  in  linear  topology  and  served  by  the  APOLLO  DN  10,000  workstation 
as  a  host  machine  (Figure  1  of  [33]).  For  the  2.62-cm  (nominal  1-in.)  model,  the  parallelized  ver¬ 
sion  of  the  code  is  about  six  times  faster  than  the  sequential  version  using  APOLLO  10,000.  It 
will  take  a  larger  memory  than  presently  available  with  this  transputer  (42  Mbytes)  to  run  a  1.31- 
cm  (nominal  1/2-in.)  model  of  the  human  body. 


WORK  PROPOSED  FOR  NEXT  YEAR 

A.  With  the  high-resolution  procedures  that  have  now  been  developed,  it  has  become  imperative 
that  we  develop  a  higher-resolution  model  of  the  human  body  than  the  one  that  we  are 
presently  using.  It  may  be  recalled  that  the  present  model  was  developed  5-6  years  ago  from 
single-tissue  properties  prescribed  for  each  1.31  cm  x  1.31  cm  area  of  the  anatomic  cross 
sections  that  were  available  for  separations  on  the  order  of  2. 3-2.7  cm.  Interpolations  were 
then  used  to  obtain  the  compositions  of  the  numerical  cells  for  various  intermediate  locations. 

With  the  possibility  of  using  linearly  expanding  grids,  we  would  like  to  be  able  to  pre¬ 
scribe  properties  for  cell  sizes  as  small  as  1.44  mm  for  dosimetric  calculations  at  3,000  MHz. 

We  would  like  to  use  a  magnetic  resonance  imaging  (MRI)  tape  of  the  human  body  to 
define  the  contours  of  the  various  organs  and  tissues  for  each  of  the  cross  sections  of  the 
body.  If  the  model  were  kept  flexible  enough  so  that  a  dominant  tissue  could  be  prescribed 
for  each  of  the  cells  likely  to  be  on  the  order  of  a  millimeter  or  so  at  higher  frequencies,  we 
could  use  a  look-up  table  of  only  16  tissues  for  all  of  the  cells  rather  than  the  presently  used 
compositions  of  the  various  tissues  for  each  of  the  N  cells  representing  the  body.  In  addition 
to  an  accurate  representation  of  the  body  at  higher  frequencies  where  cell  sizes  on  the  order  of 
millimeters  are  needed,  use  of  a  16-tissue  look-up  table  would  result  in  a  considerable  reduc¬ 
tion  in  the  memory  requirement  and  hence  the  computation  time  for  the  code. 

B .  We  would  like  to  investigate  techniques  to  improve  efficiency  of  calculations,  particularly  at 
higher  frequencies  where  the  cell  size  is  of  necessity  small  and  hence  a  large  number  of  cells 
are  needed  to  represent  the  whole  body. 

At  higher  frequencies  where  the  EM  deposition  is  primarily  superficial,  it  may  be  pos¬ 
sible  to  bring  the  absorbing  boundary  behind  the  model  closer  to,  and  perhaps  within,  the 
body,  thereby  reducing  the  interaction  space  that  is  to  be  modeled. 
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A  method  of  choice  will  be  the  use  of  the  linearly  expanding  grid  by  which  the 
exposed  frontal  region  of  the  body  may  be  expressed  by  a  finer  grid  and  an  increasingly 
coarser  grid  used  for  the  interior  of  the  body.  Now  that  this  algorithm  has  been  developed 
and  tested  with  objects  such  as  spheres  (Figure  5),  it  should  be  possible  to  adapt  it  for  use 
with  the  new  high-resolution  model  of  the  human  body  that  will  be  developed  as  Task  A  of 
the  project. 

Once  the  high-resolution  anatomically  based  model  of  the  human  body  has  been  devel¬ 
oped,  it  will  be  our  objective  to  develop  an  atlas  of  SAR  distributions  for  frequencies  up  to  5- 
6  GHz  using  the  techniques  outlined  here.  It  would  thus  be  possible  to  obtain  SAR  distribu¬ 
tions  for  the  biologically  sensitive  regions  such  as  the  eyes,  the  brain,  the  gonads,  and  oth¬ 
ers,  w'ith  resolutions  on  the  order  of  millimeters. 

For  this  atlas  of  SAR  distributions  for  isolated  and  for  grounded  and  shoe-wearing 
conditions,  we  would  cover  the  region  10-6,000  MHz  with  smaller  frequency  separations  of 
10  MHz  at  lower  frequencies  increasing  to  as  much  as  200-500  MHz  at  the  high  end  of  the 
frequency  band.  We  have  talked  to  the  personnel  at  the  Utah  Supercomputing  Institute 
(University  of  Utah)  who  have  volunteered  to  help  us  develop  the  computer  graphical  proce¬ 
dures  to  display  the  calculated  SAR  distributions.  When  these  procedures  are  developed  it 
should  be  possible  to  include  computer-generated  pictures  of  SAR  distributions. 

C.  Once  Tasks  A  and  B  have  been  satisfactorily  completed,  we  would  like  to  use  the  technique 
thus  developed  to  understand  the  coupling  of  the  high-power  microwave  (HPM)  pulse  to  the 
human  body.  The  exposure  parameters  to  use  for  these  calculations  (frequency,  time 
domain,  wave  shape,  etc.)  will  be  prescribed  by  Dr.  Erwin  of  the  USAF  School  of 
Aerospace  Medicine. 

If  pertinent,  the  nonlinearities  of  the  tissue  properties  at  the  high-field  strengths  associ¬ 
ated  with  HPM  pulse  may  also  be  incorporated  by  modifying  the  Maxwell's  equations  that 
are  solved. 

D.  The  data  calculated  for  coupling  of  the  HPM  pulse  will  be  graphically  displayed  for  the  vari¬ 
ous  sections  of  the  body  by  means  of  a  videotape  of  a  computer-generated  movie. 

E.  The  impedance  method  previously  developed  in  our  laboratory  [27-29]  is  ideally  suited  for 
calculations  of  induced  currents  and  SARs  for  exposure  to  time-varying  magnetic  and  electric 
fields  due  to  electromagnetic-rail-gun  type  devices.  We  look  forward  to  receiving  the  specifi¬ 
cations  of  such  a  device  from  the  Air  Force  project  monitor  to  proceed  w  ith  the  calculations 
for  the  high-resolution  anatomically  based  model  of  the  human  body. 

F.  Significant  progress  has  recently  been  made  in  developing  an  anatomically  based  thermal 
model  of  the  human  body  [ 30)  to  calculate  temperature  distributions  for  EM  exposure  condi¬ 
tions.  In  this  work,  a  computerized  tomography  (CT)  scan-based  model  of  the  human  torso 
has  been  combined  with  the  anatomically  based  model  for  the  rest  of  the  body  to  obtain  tem¬ 
perature  distributions  in  the  presence  of  SARs  due  to  RFR.  The  model  of  the  body  consisted 
of  about  150,000  cubic  cells  (0.87  cm  per  side),  and  about  450,000  air  cells.  Unlike  the 
previous  model,  this  model  incorporates  a  skin  layer  which  allows  for  a  proper  heat  balance 
with  the  surrounding  ambient.  We  propose  to  adapt  this  model  to  the  SAR  calculations  that 
are  now  possible  to  about  3,000  MHz  and  start  to  obtain  the  temperature  distributions  for  the 
RFR  exposure  conditions  at  various  frequencies.  Of  particular  interest  are  the  power  densi¬ 
ties  (or  SARs)  that  will  result  in  measurable  changes  in  temperature  for  the  various  organs 
and  the  ensuing  changes  in  the  biological  parameters  such  as  blood  flow  and  heat  loss  by 
mechanisms  such  as  perspiration,  respiration,  and  others. 
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CONCLUDING  REMARKS 


The  major  objective  of  this  project  is  to  develop  advanced  user-friendly  codes  with  degrees  of 
resolution  that  would  have  been  unthinkable  just  a  few  years  ago.  The  eventual  goals  of  the  project 
are  to  create  an  atlas  of  SAR  distributions  for  the  human  body  for  RFR  exposure  for  isolated  and 
for  grounded  and  shoe-wearing  conditions,  and  the  fifth  edition  of  the  RFR  dosimetry  handbook 
with  emphasis  on  anatomically  based  models  for  SAR  and  temperature  distributions. 
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